Controlled Flow of Spin-Entangled Electrons via Adiabatic Quantum Pumping 
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We propose a method to dynamically generate and control the flow of spin-entangled electrons, 
each belonging to a spin-singlet, by means of adiabatic quantum pumping. The pumping cycle 
functions by periodic time variation of localized two-body interactions. We develop a generalized 
approach to adiabatic quantum pumping as traditional methods based on scattering matrix in one 
dimension cannot be applied here. We specifically compute the flow of spin-entangled electrons 
within a Hubbard-like model of quantum dots, and discuss possible implementations and identify 
parameters that can be used to control the singlet flow. 
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Entanglement is one of the most intriguing features 
that distinguish the quantum world from the classical. 
In recent years Q applications potential in quantum in- 
formation and computation has injected renewed vigor 
into the study of entangled states. Since Bohm's [2| refor- 
mulation of the Einstein-Podolsky- Rosen (EPR) paradox 
Q in terms of pairing of spins, spin singlets have become 
the canonical example of an elementary entangled state. 
From an experimental standpoint, manipulating singlet 
pairs of electron spins, in particular, is promising because 
of the vast expertise already available in solid-state elec- 
tronics 

Several proposals have emerg ed HH00 that aim to 
generate a controlled flow of electron spin singlets using 
Coulomb blockade and tunnelling at pinched quantum 
dots that have weak coupling to leads. In these pro- 
posals, the natural tunnelling of unwanted single elec- 
trons can be difficult to suppress and/or singlet delivery 
time is constrained by tunnelling rates that are slow due 
to the necessity for pinched dots. Since manipulations 
in a quantum information device need to occur before 
decoherence sets in, fast delivery is essential, and any 
background of unpaired electrons reduces the purity of 
entanglement. It is therefore highly desirable to have 
a fast dynamic mechanism to generate a selective flow 
only of electrons belonging to singlets. In this paper we 
show that this can be accomplished through a generalized 
adiabatic quantum pumping that is induced by localized 
interactions and we discuss possible ways of experimen- 
tal realization. Pinched dot tunnelling times and singlet 
formation times do not constrain delivery rates in our 
proposal. In addition, the absence of bias and the adia- 
batic nature of the pumping can reduce the production 
of heat in the apparatus. 

The notion of quantum pumping has its roots in a spec- 
ulative paper by Thouless [9J in 1983, but advances in 
nanoscale transport have led to a renewed and growing 
interest in the phenomenon in recent years both theoret- 
ically and experimentally El El El El El El E3 
Quantum pumping is a coherent process that creates 
a direct current in the absence of any bias through a 



nanoscale device, by changing its scattering properties 
periodically through independent adiabatic variation of 
two or more physical parameters. Adiabatic quantum 
pumping of charge E3; s P m Ql an d thermal E3 
currents have been considered. However previous studies 
have generally relied on a theoretical description based on 
transmission and reflection coefficients which cannot be 
used to describe the pumping of singlets due to localized 
interactions; therefore in this paper we also formulate 
a more generalized approach to the theory of quantum 
pumping. 

Singlet Current: We consider a single available chan- 
nel in a quasi one dimensional mesoscopic conductor con- 
nected to macroscopic contacts (Fig. da)). In the pres- 
ence of two body interactions the charge current can be 
defined in terms of the two particle reduced density ma- 
trix, p 2 . If the interaction does not affect spins, the 
reduced density matrix separates into four independent 
spin subspaces, one singlet with a symmetric spatial part 
and three triplets with antisymmetric spatial parts; the 
singlet current is therefore equivalent to the charge cur- 
rent associated with the symmetric spatial part (pf ) of 



V(xi, x 2 ; t) 
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FIG. 1: (a) Schematic figure showing the iwo-body inter- 
action V{xi,X2\t) acting in a finite interval (—1,1) in a one 
dimensional system, (b) Implementation in a tight-binding 
model where the two-body interaction is present only at lat- 
tice sites m = ±1; the interaction strengths U±(t) at the two 
sites are the time-dependent pumping parameters. 
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the two particle reduced density matrix 
eh f 

J s {xi,t) = — / dx 2 Q{d Xl p 2 (xi,x 2 ;x' 1 ,x 2 ;t)} x , .(1) 

171 J l 

We consider a situation in which two-body interactions 
are localized, meaning that they are non-vanishing only 



when the particles are in a certain finite interval, as 
shown in Fig. |TJa). Such interactions lend themselves 
to a scattering description and the current can be ap- 
proximately evaluated by expanding the density matrix 
in terms of two particle scattering states 



J s ( Xl: t)^^- JdEF(E)J ^/^f + ^Sf ~^)/ fe^i^^^,^^,^,*)^^,^^,^^)}^^^ ■ (2) 

I 



Here E denotes the energy required to remove a pair of 
particles from the many-body ground state, and F(E) 
is the distribution of this pair energy. The effect of in- 
teraction on the current is determined completely by the 
two-particle singlet scattering states, ^k 1: k 2 arising from 
free singlet states $ kltk2 (xi, x 2 ) = -^[<t'k 1 (xi)(t>k 2 (x 2 ) + 

0fci {x 2 )4>k 2 {xi)] where 4>k{x) denotes a single particle 
plane wave state with momentum hk. 

Pumped Current through Adiabatic Perturbation : The 
two particle scattering states, and therefore the current, 
are determined by the interaction V(x, t) between a pair 
of particles; we take it to be time-dependent and to occur 
only in a finite region \xi\ < I. Most importantly the in- 
teraction V(x,t) is chosen to be localized so that it only 
affects singlets thereby naturally eliminating the flow of 
triplets in the absence of a bias. When the characteristic 
period ix> of the time variation of the potential is slow 
compared to the time St the particles dwell in the scat- 
tering region [17j, w x St <C f, we can apply adiabatic 
perturbation theory to express the scattering states of 
the time-dependent Hamiltonian in terms of the instan- 
taneous states up to linear order 

*- k (x,t)~mx)-ih J dx'G* (x, x'; E) . (3) 

We use the notation x = {xi,x 2 } and k = {ki,k 2 }, 
so that G t {x,x';E) is the iwo-particle instantaneous re- 
tarded Green's function for the full Hamiltonian. The 
instantaneous state ty%(x) is a solution of the time- 
independent Lippmann-Schwinger equation for the po- 
tential V(x',t) at the specific time t, 

= * s (s) + J dx'G (x,x';E)V(x',t)^- k (x')(A) 

where Gq is the free two-particle retarded Green's func- 



tion. Taking the time derivatives of the defining equa- 
tions for \&*(af') and G*(x,a;';-E) enables us to express 
the second term in Eq.(|3J|, that is linear in dt, as 

A* £ (x,t) = -ih J dx J dx"G t (x,x';E) (5) 

xG t (x',x";E)V(x",t)mx"). 
If there is no bias or time-dependence, the laws of 
thermodynamics demand that there should be no cur- 
rent; we explicitly confirm that our expression for the 
current satisfies this essential physical requirement. The 
net current in the absence of time-dependence is eval- 
uated by using the zeroeth order term from Eq. @ 
for the scattering state in Eq. J2J: ^k 1 .k 2 ( x i, x 2 ,t) = 
^k k ( x i,x 2 ). The resulting expression can be simplified 
by relating the imaginary part of the retarded Green's 
function to the free singlet states: 9{Go(x, x'; E)} = 

-nJdhJdk.S (™L + ™L-e) <S> k (x)<Pl(x>). Then 
repeated use of the Lippmann-Schwinger equation and 
properties of the Green's functions shows that the net 
current corresponding to the zeroeth order 4 ,t (x) van- 
ishes. 

After confirming that our expression cannot produce 
spontaneous current, we evaluate the singlet current in- 
duced by the adiabatic time evolution. To linear order 
in the time dependence this involves the evaluation of 

J dx^id^Kxux^Amx^x^t)}^^ (6) 

within the expression for the current in Eq. p|. A cal- 
culation employing standard Green's function identities, 
similar to that for the zeroeth order, leads to an expres- 
sion for the net amount of singlet entangled electron pairs 
pumped in a complete cycle of period r, yield the main 
result of this paper: 



Qs{ r) = ^-[dtJdEF { E)^[Jd^,t)J 



z 2 \s 



{G t *{x,x';E)d Xl G t (x,x';E)} 



(7) 
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Singlet Pumping in a Turnstile Model: We illustrate 
our results with a tight-binding model, (Fig. [Jb)) with 
two Hubbard impurities located at sites —to, m 

(8) 

An electron can interact with another only at those two 
sites therefore, due to the Pauli principle, only singlets 
are affected. The strength of the interactions, U±(i) are 
the two time-dependent pumping parameters. This con- 
cept is similar to a 'turnstile model' but differs sig- 



nificantly in that, instead of time- varying external po- 
tentials, the two-body interaction among electrons is var- 
ied in time. We separately derived a discrete version of 
Eq. JJJ); the end result amounts to replacing the coordi- 
nate arguments with site indices x — > n, integrals by sums 
and derivatives with a finite difference form. A lengthy 
calculation leads to an expression for the singlets pumped 
in a complete cycle in terms of the free two-particle lattice 
Green's function, specifically two of its matrix elements 
G (0) = G (m, to; E) and G (2m) = G (m, -to; E) 



, -e / d STr-T M |T± W|2 NGoWH 1 + \T T (t)G (2m)\ 2 ) ± 25{r T ft)G (2m)G Q ± (2m)}] 
!s(t) = — I dtjdKb{h)—^U±(t) TTTTvIu t nw u\r< fo^\\2 • ( 9 ) 



2tt 



U ± (tn\ - T T (t)T ± (t)G (2m)G Q (2mW 



r 



Here m = {m,m}, and T±(t) = l/(U ± \t) + G o (0)) 
is the T-matrix for a single Hubbard impurity, and 
G^ = Go (2m) for A = ±i|A|, when expressed in the form 
Go (2m) = /; f with cosh(A) = E/2 - cos(fc). 

Exact analytical forms exist for the lattice Green's 
functions, Go, in terms of elliptic integrals ^^|. For 
the purpose of numerical estimates, we assume a square- 
profile time dependence in the plane of the parame- 
ters U±(t), shown in Fig.|2fa), where the two parameters 
change alternately between minimum value V m m and a 
maximum value U ma x- The pair distribution function is 
taken to be a Fermi function, F(E) ~ f/[e /3(B ~ £) + 1]. 
At low temperatures, an integration by parts with respect 
to energy yields an expression for the pumped singlets in 
terms of the maximum energy, £ available for a pair. 

The onsite energy of each tight-binding site is taken 
to be zero and all energies are expressed in units of the 

(a) 



Qs(T) 




2 U max 



FIG. 2: (a) Pumping cycle in the space of parameters U±. 
(b) Net singlets pumped in a single cycle as a function of the 
location and size of the square pumping cycle as Umin and 
Umax are varied. 



nearest neighbor coupling strength. In Fig. |2Ib) we plot 
the net singlets pumped in a single cycle as a function of 
the size and location of the square footprint of the time- 
cycle in the space of the parameters U±; the flow depends 
on the enclosed region. Figure |3| shows the dependence 
of the singlet current on the parameter £ that measures 
the available energy for pairs of electrons determined by 
the chemical potential in the contacts. The two curves in 
the figure correspond to different locations of the Hub- 
bard impurities, at lattice sites m = ±1 and m = ±2, 
illustrating the significant effect the spatial separation of 
two impurities have on the pumping rate. The direction 
of flow can also reverse for certain values of the various 
parameters; reversing the time-cycle is not the only way 
to reverse the direction of the current [loj . Quantum 
pumping has the intrinsic property that the magnitude 
of the pumped quantity, in this case singlets, is continu- 
ous in nature so that the delivery rate per cycle can be 
continuously adjusted. Thus there are several ways to 
precisely control the magnitude and direction of the flow 
of singlets dynamically. 

Discussion and Outlook: The turnstile model could 
be implemented by taking the interaction sites to be 
quantum dots coupled to leads and varying the electron- 
electron interaction by decreasing the size of the dot pe- 
riodically. Concurrent variation in the dot-lead coupling 
could in principle be compensated by counter-varying 
gate voltages. Another way to vary the interaction 
strength is to change the dielectric constant locally. This 
could be accomplished by a local shift of the electron 
wavefunction among layers of different materials in a het- 
erostructure; this technique was recently used to control 
the electron g factor [Pj. A shift between silicon and 
germanium layers or gallium arsenide (GaAs) and alu- 
minium arsenide (AlAs) layers can in principle produce 
up to 25% change in the dielectric constant, sufficient to 
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FIG. 3: Net singlets pumped per cycle as a function of the 
maximum available energy £ for pairs. The two curves cor- 
respond to different separations of the Hubbard impurities, 
solid line — > ±m = ±1 and dashed line — > ±m = ±2. 



see the effects described here. 

The adiabatic condition requires the period, r 3> St, 
the dwell time given by St = d/v, with d being the size of 
the scattering region and v the carrier velocity ^3 • Since 
our singlet pump requires no potential barriers that re- 
duce kinetic energy, v may be taken to be the Fermi veloc- 
ity typically > 10 5 m/s, so for a scattering region d ~ 10 
nanometers the adiabatic condition would allow thou- 
sands of cycles per nanosecond. Our numerical estimates 
then give a pumping rate at the order of a thousand sin- 
glets per nanosecond. In typical Coulomb-blockade based 
schemes the most optimistic estimates yield a delivery 
rate at the order of one nanosecond per singlet 0, 0- 
Thus our approach has the potential to be much faster. 

Our result, Eq. Q, has the merit that it can also be 
applied to quantum pumping in systems that allow an in- 
dependent particle description, as considered in previous 
studies. All the elements in Eq. Q then reduce to single 
particle functions: F(E) — > f(E) is the Fermi distribu- 
tion function, the Green's function is a single particle 
one with an asymptotic form — > — i[m/ \Ti 2 k)]e lkx i{) t t ?(x') 
in terms of the scattering state ipl.(x') and wavevector 
k = [2mE/h ]~ 1 ^ 2 . This yields the pumped charge 

which agrees with expressions derived in earlier works 
But unlike all previous treatments, we never use 
ID scattering matrix elements as they do not have useful 
generalizations when particles interact with each other in 
a region rather than scatter off an external potential. 

To summarize, we have proposed a method based on 
two-body adiabatic quantum pumping for generating a 
dynamically controlled flow of spin-entangled electrons. 
The process is inherently coherent, potentially much 
faster than most current proposals, has reduced noise 
because of the lack of bias and the natural elimination of 
single electron flow, and allows for continuous adjustment 
of the flow through numerous physical parameters. All 
of these features can be developed and incorporated into 



a comprehensive scheme to generate a controlled flow of 
entangled electrons. Our goal here has been to present 
the basic idea, develop a theoretical framework for its 
description and discuss a possible physical model for im- 
plementation. We have in the process generalized the 
treatment of quantum pumping to incorporate interac- 
tions that cannot be treated in a scattering matrix ap- 
proach. We anticipate that future studies can build upon 
the considerations in this work. 

This work was supported by the Packard Foundation 
and NSF NIRT program grant DMR-0103068. 

Note added: When we were finishing this paper we be- 
come aware of an excellent recent proposal by Beenakker 
et al. iii to excite entangled electron-hole pairs using 
one-body potential in a scattering matrix approach. Our 
proposal is more difficult to implement, but it could be 
naturally incorporated into an electron-only spintronic 
device, it produces a current of singlets unadulterated by 
single electrons, and it will enjoy a longer decoherence 
time since hole decoherence times tend to be short. 
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